Spontaneous and resonant lifting of the spin blockade in nanowire quantum dots 
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A complete description of the dynamics of a two-electron system confined in narrow nanowire 
quantum dots under oscillating electric field and phonon mediated relaxation is presented in the 
context of recent electric dipole spin resonance experiments. We find that spin-orbit coupling results 
in lifting the spin blockade by spin non-conserving phonon mediated relaxation provided that the 
initial state is close in energy to the ground state which results in deblocking one of the Zeeman split 
triplets. At higher magnetic field, after singlet-triplet anticrossing new channel for lifting the Pauli 
blockade opens which results in an appearance of additional resonance lines, as present in recent 
experimental results [S. M. Frolov, et al., Phys. Rev. Lett. 109, 236805 (2012)]. 
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In recent years gate defined nanowire [l[ double quan- 
tum dots have been successfully used for experimental 
demonstration of electrical control of single spins 

0| . The spin rotations are performed by means of elec- 
tric dipole s pin resonance (EDSR) where spin-orbit (SO) 
interaction [lC 



ll| is used to electrically control the spin 



excluding the need for oscillating magnetic field in the 
device [12[. The spin oscillations are probed exploiting 
the spin blockade [13] of two-electron system where the 
current cycle (0,1) (1,1) (0,2) (0,1) [the num- 
bers denote number of electrons in the adjacent quantum 
dots] is blocked at the (1,1) — » (0, 2) transition when the 
spin configurations of the (1,1) and (0,2) states do not 
match. For low bias the only available (0,2) state is 
the spin singlet so the current is blocked if the system 
is initialized in one of the spin-polarized triplets. The 
blockade is lifted when the spin of oscillating electron is 
rotated such the total spin of the (1,1) state is changed 
from S = 1 to S = 0. Relaxation of the (1,1) state 
that follows the spin rotation opens the tunneling of the 
single electron through the double dot and is mediated 
by phonons that disperse the excess energy between the 
(1,1) and (0,2) states. 

Strong SO coupling in InSb and In As nanowires that 
is useful for effective control of the spins leads to spin re- 
laxation 14-16[ that results in lifting the spin blockade, 
limits the available magnetic fields for the EDSR [I?} and 
results in the leakage current that exhibit dips or peaks 
at low magnetic field depending on the interdot coupling 
[l8| . In the present Letter we study phonon mediated 
electron relaxation in electrically driven nanowire quan- 
tum dots. We find several unexpected features that are 
crucial for the mechanism of the spin blockade lifting: i) 
in all regions where the spin conserving (1,1) — >• (0, 2) re- 
laxation occurs the spin non-conserving relaxation from 
(1,1) triplet state with spins polarized along the mag- 
netic field is of a very similar effectiveness, ii) For small 
magnetic fields where (0,2) singlet is the ground state 
this leads to the spontaneous lifting of the spin blockade 



from this triplet. Hi) On the other hand the relaxation 
from triplet with spins oriented antiparallel to the mag- 
netic field orientation is two orders of magnitude slower 
so the blockade is maintained, iv) At higher magnetic 
fields when the triplet becomes the ground state spin ro- 
tation accompanied by charge redistribution results in 
lifting spin blockade through direct transition to (0,2) 
singlet which is visible in recent experimental maps [5j. 

The considered two-electron system is described by 
the Hamiltonian H(t) = ^^(t) + e 2 /(47ns£ |ri - r 2 |) 
where h l (t) is single electron energy operator. For a 
narrow nanowire the charge occupies ground state of 
lateral quantization (here taken in the Gaussian form 



i/>(y,z) = (y/jTL 
which leads 



da tc 



exp[-(?/ 2 + z 2 )/2l% with I = 20 nm) 
to the two-electron Hamiltonian, 



H lD (t) = h\ D {t) + h\ D {t) + ^^erf 
with single-electron energy operator 



\xi - x 2 \ 



a/2/ 



(1) 



h 1D (t) 



acF y k x 



-H B g(x)B(j x , (2) 



where hk x = —ihV x is momentum operator and Hso — 
—aciykx stands for Rashba SO coupling. We allow for 
position dependent ^-factor in the device [H, 0, H| and 
take g(x) = g[l + (3H(x)] where H(x) is Heavyside step 
function and f3 = 0.1. The potential is separated into 
V(x,t) = Vqd(x) + V'(x,t) where quantum dot confine- 
ment is described by Vqd(x) = V c (x) + eFb ias x. V c de- 
fines potential of two quantum dots of 138 nm width each 
separated by potential barrier of 25 nm width and 200 
meV height. The driving AC electric field is assumed ac- 
tive in the left dot [2[ , so the time dependent part of the 
potential takes the form V f (x, t) = eFAcxf(x) s'm(u Act) 
where f(x) = 1 in the left dot and outside - see the 
inset to Fig. QJa). 

We calculate eigenstates ^ n (xi, ai, x^, c>"2, t = 0) 
(with corresponding eigenenergies E n ) of Hamilto- 
nian Eq. (pQ) for t = and use them as 
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FIG. 1. (color online) (a) Energy levels of two-electron double 
quantum dot in function of bias electric field for B — 50 mT. 
The arrows presents approximate spin polarization of elec- 
trons in the dots. The inset present schematics of the consid- 
ered confinement potential, (b) Half-life time of excited states 
due to phonon mediated relaxation to the (0,2) singlet state. 
The colors of curves denote initial state of relaxation. The 
symbols (curves) corresponds to the results obtained without 
(with) SO interaction. 



a basis for time evolution [20| where the two- 
electron spinor is expressed as ^(xi, cr l5 X2, <?2, t) = 
En °n(t) exp(-iE n t/h)^ n (x u a u x 2 ,(J2,t = 0). 

To describe the phonon mediated relaxation we allow 
for the transitions between the two-electron states with 
rate given by the Fermi golden rule. The relaxation rate 
between the initial and final states is described by, 



--1 = _ V 

if h 

x |(^|e 



dq\M„(q)\ 2 
^\^)\ 2 5{\E f 



(3) 



where the phonon dispersion relation is E q = ?kv|q| and 
c v is the sound velocity. The sum in Eq . (Ip goes over 
three types of electron-phonon scattering [20| due to: de- 
formation potential with longitudinal mode [y =LA-DP) 
(2lj . piezoelectric field (y =LA-PZ) and piezoelectric field 
with transverse modes (y =TA-PZ) [22j. 

The relaxation is included into time dependent calcu- 
lation such the amplitude \c n (t)\ 2 of each eigenstate is 
changed due to relaxation to all lower energy states and 
is increased due to relaxation from higher energy states 
with corresponding transition rates r^j [io|. This cor- 
responds to OK temperature where the energy instantly 
disperses and results in a half-life time of an i'th state de- 
fined as T1/2 = \n(2) /rif due to relaxation to /'th state. 

We assume material parameters for InSb, i.e. electron 
effective mass m* = 0.014, g = —51, dielectric constant 



e = 16.5 and take the Rashba constant a = 10 meVnm. 
The AC field amplitude Fac — 0.05 kV/cm is assumed. 

The charge distribution in the double dot is controlled 
by external voltages applied along the structure. Fig. 
[T](a) presents the lowest part of the energy spectrum (the 
subsequent energy levels - of (0,2) triplets - are above 5 
meV) in function of bias electric field for B = 50 mT. For 
the most negative values of Fbi as the ground state is a sin- 
glet |*,|t — ti) state for which both electrons reside in 
the right dot [(0,2) configuration]. The four excited states 
correspond to single occupancy of each dot [(1,1) config- 
uration] and energy of those states only weakly change 
in function of bias electric field. The two close in en- 
ergy states, i.e. | |, t) and | t, |) correspond to definite 
and opposite spin configurations in each dot (i.e., in the 
I I, t) state spin of the electron in the right dot - where 
the ^-factor takes the highest value - is polarized along 
the magnetic field) resulting from mixing spin-zero triplet 
with singlet state by the ^-factor mismatch between the 
dots. The two triplets | t>t)j I are s P n ^ be Zeeman 
interaction. 

Let us inspect relaxation times of the four excited (1,1) 
states. In Fig. [U^b) we present half-life times of excited 
states due to relaxation to the ground state singlet. In 
the absence of SO interaction phonon scattering couples 
only states with the same total spin. Only the half-life 
times of | t, |) and | |, t) have finite values and they are 
presented with the crosses in Fig. [ljb). At low values 
of Fbi as the times are of order of hundreds of microsec- 
onds but when the energy differences between the initial 
states and (0,2) singlet become lower the half-life times 
rapidly drop allowing for (1,1) — » (0,2) spin-conserving 
relaxation within nanoseconds for Fbi as > —0.16 kV/cm. 

When SO coupling is included the spin polarization 
of the states becomes only approximate. The relaxation 
times of | |, t) and | t, |) states do not change - see 
the curves and crosses in Fig. QJb). However now re- 
laxation from all the (1,1) states to (0,2) singlet is open. 
For most negative values of Fn as half-life time of triplet 
states is longer than millisecond but when the bias field is 
increased the half-life time of | t, t) triplet becomes about 
the same as the two spin opposite states lifting the spin 
blockade. Only the | l\) state remains effectively blocked 
with relaxation time longer by two orders of magnitude 
from the rest of the (1,1) states. This scenario was al- 
ways repeated as we checked for different strengths of SO 
coupling and the dot size [2o| . 

The experimental studies [7] report spin coherence 
time of order of tenths nanoseconds and coherent manip- 
ulation over a single spin up to 100 ns. We therefore focus 
on Ffcas range where the relaxation times are of order of 
tenths of nanoseconds - hereafter we take Fbi as = —0.15 
kV/cm - that allow for deblocking of single-electron cur- 
rent through the double dot after spin rotation in EDSR 
experiments. In Fig. EJa) we plot energy levels in func- 
tion of the magnetic field B. For B = the excited state 
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FIG. 2. (color online) (a) Energy spectrum in function of 
the magnetic field. Straight arrows denote transitions due 
to phonon relaxation. Curved arrows depict available EDSR 
resonances from the triplets, (b) Half-life time of excited 
states due to phonon mediated relaxation. Results obtained 
for F bias = -0.15 kV/cm. 



is fourfold degenerate due to high interdot barrier - neg- 
ligible exchange coupling. When the magnetic field is in- 
creased the energy levels of the two spin polarized triplets 
| t> t) and | i, 1} are split by the Zeeman interaction. On 
the other hand the energy levels of the two spin-opposite 
states - | |, t) and | t, 1} ~ are weakly split due to ^-factor 
mismatch in the dots. At B = 0.1 T an anticrossing be- 
tween the energy levels of triplet | t) and (0, 2) singlet 
states appear followed by the change of the ground state. 

Figure [2fb) presents half-life times of excited states 
due to relaxation to (0,2) singlet. We observe that the 
relaxation from spin antiparallel | ^, \) and | 1} states 
occurs within few nanoseconds regardless of B value. As 
the energy separation between energy levels of | t, t) and 
(0,2) singlet decreases the relaxation time drops and af- 
ter B = 50 mT the half-life time of this state is even 
lower than the half-life time of spin-antiparallel states. 
On the other hand relaxation from | ^, \) state is slow 
and the half-life time grows for increasing magnetic field 
until B = 0.1 T. This shows that for magnetic field range 
before the anticrossing only the | |,|) triplet provides 
spin blockade as the | t?t)> I -Jot) and | states decay 
quickly into |»,|t - t|). 

In EDSR experiments the two-electron system can ini- 
tialize in any of the low energy states within the transport 
energy window. We therefore study the time evolution 
taking each of the (1,1) states with the energy levels de- 
picted in Fig. [2fa) as an initial state. In Fig. [3je)-(h) we 
present (0,2) occupation probability (which would allow 



for tunneling of one of the electrons outside the dot lift- 
ing the blockade) averaged during 30 ns time evolution 
in function of the magnetic field and the electric field 
frequency. 

For | tt) taken as the initial state the probability is 
presented in Fig. [3fe). At the left part of the map we 
observe increasing probability in the background in func- 
tion of B due to spin relaxation that results in sponta- 
neous lifting of spin blockade. At B = 0.1 T the phonon 
relaxation to singlet (0,2) from | t, t) stops as the latter 
becomes the ground state - the background of the plot 
shows nearly zero (0,2) occupation probability. However 
we observe several resonance lines with increased prob- 
ability. The resonance line (A) corresponds to the spin 
rotation in the left dot (| t>t) I I A)) accompanied 
by the phonon mediated relaxation to the |»,|t — tl) 
singlet which results in an increase of the (0, 2) occupa- 
tion probability - see Fig. [3]^a) where we present the 
probability \c n \ 2 of finding the system in the n'th state 
during the time evolution. The (0) transition is related 
to the spin rotation in the right dot which is much less 
effective due to presence of the AC electric field only in 
the left dot and high interdot barrier which results in a 
narrow resonance line. The bottom line marked with (□) 
corresponds to the direct transition to the (0, 2) singlet 
that involves charge reconfiguration between the dots - 
see Fig. [3^b). Note that line of increased probability due 
to | t>t) ~~ ^ | # ?it — tl) transition is not observed for 
B < 0.1 T as in this region the spin relaxation of the 
triplet results in its fast deexcitation to the ground state 
with rate that exceeds the EDSR transition. 

In Figs. [3jf), (g) one can observe mainly nonzero (0,2) 
occupation probabilities due to fast spin-conserving re- 
laxation of | |,t) and | t>l) to the |»,|t - ti) state 
as discussed previously. In fact this relaxation is fast 
enough that one can observe lines of lowered probability 
when the system already relaxes into singlet (0,2) and is 
driven back to one of the excited states. 

For the | |, |) triplet taken as the initial state outside 
the resonances the (0,2) occupation probability is nearly 
zero at Fig. [3fh) as the phonon mediated relaxation from 
this state is slow - see the blue curve in Fig. [2jb). The 
lines that go through the diagonal of the plot - (♦), (A) 

- corresponds to the transition to the | |, t) and | t, \) 
states respectively accompanied by relaxation to |*,|t 

— ti) [see Fig. [3f d)] and the line at the left upper part 
of the plot - (■) - is a direct transition to the (0, 2) 
singlet that does not involve phonon mediated relaxation 
[see Fig. |3(c)]. 

Note that in maps of Fig. [3ja) and (d) also lines of 
increased probability at the half frequency of the (□) 
and (■) transitions are visible which is due to resonant 
harmonic generation by the driven electrons [23|. 

The experimentally [5[ measured resonances at current 
maps are obtained from many sequential events of single 
electron transport through the structure. In each of them 
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FIG. 3. (color online) (a)-(d) transitions between the eigenstates during the time evolution at the resonances marked with the 
symbols, (e)-(h) Probability of (0,2) occupation averaged during the 30 ns time evolution obtained for subsequent (1,1) states 
taken as the initial state of the time evolution - the spin configuration of the initial state is denoted in the left bottom corner 
of each plot. 




FIG. 4. (color online) Probability of (0,2) occupation aver- 
aged over the 30 ns time evolution calculated as a sum of 
results obtained for initial states with (1,1) occupation. For 
each B value the (0,2) occupation probability obtained for 
huoAC = (in the absence of driving electric field) was sub- 
tracted from the results. 



the system can initialize in any of the (1,1) states. We 
therefore calculate the total probability of (0,2) occupa- 
tion by summing the results for initial states presented 
in Figs. Efe)-(h). For each value of B the probability 
obtained without (due to pure relaxation) the oscillat- 
ing electric field is subtracted to mimic the experimental 
procedure of Ref. [5] . For low values of B we observe two 
lines at the diagonal of the map that corresponds to the 
transitions from | , \) state - rotation of the spin down 
to spin up in the left dot (bright line) or in the right 
dot (faint line) accompanied by relaxation to (0,2) sin- 
glet. After singlet-triplet anticrossing at B = 0.1 T the 
lines correspond to transition from both the triplets. At 
B = 0.1 T additional resonance line starts at the bottom 



of the plot that corresponds to spin rotation with charge 
reconfiguration from | t,t) triplet. Note that there is no 
similar line corresponding to transition from the | 
state as it is compensated by the lowered probability ob- 
tained for evolution starting from | ^, t) and | t>i)- 

Similar arrangement of the resonance lines is preset in 
experimental map of Ref. (EJ that probed wider range 
on magnetic field as compared with the previous experi- 
mental studies [!-[3, 0, [?J . Especially the line that corre- 
sponds to (□) direct transition is visible which indicates 
the ground state change occurred blocking the low energy 
triplet state. One has to keep in mind however that in 
the experiment Ref. [H[ the ^-factor mismatch was com- 
pensated in some regions of the plot resulting in a single 
diagonal resonance line therein. 



In conclusion we have presented the role that spin re- 
laxation and EDSR plays on lifting the spin blockade in 
coupled nanowire quantum dots. We found that spin re- 
laxation can lead to spontaneous lifting of the spin block- 
ade such that the resonances are observed only from a sin- 
gle spin-blocked triplet state. The change of the ground 
state in higher magnetic fields leads to spin blockade of 
both the triplets and reveals additional resonance to the 
(0,2) singlet that do not involve phonon mediated relax- 
ation and which is present in recent experimental results. 
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